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The recent chirally extrapolated result of Leinweber and Thomas [Phys. Rev. D 62, 074505 (2000), 
or LT] for the nucleon strangeness form factor G^j(O) — -0.16 ± 0.18 fiN differs markedly from the 
earlier result -0.75 ± 0.30 hn obtained by Leinweber [Phys. Rev. D 53, 5115 (1996)] from the same 
lattice data. An unresolved problem in the LT analysis of lattice data is identified and addressed. 
A value of G%f{0) — -0.55 ± 0.37 /iat is obtained at i?3 = 0.55 by extrapolating only the nucleon 
isoscalar lattice data. 
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PACS number(s): 12.38.Gc,12.38.Lg,13.40.Gp 

Past lattice calculations made by Leinweber and col- 
laborators of baryon magnetic moments (MMs) have 
been restricted to massive quarks equivalent to pion 
masses m greater than 600 MeV. The calculated mo- 
ments have to be extrapolated to the physical pion mass 
= 140 MeV before comparison with experiment. Ex- 
trapolating linearly in m^, Leinweber has estimated 
the nucleon strangeness form factor at momentum trans- 
fer Q2 = to be 01/(0) = -0.75(30) in units of 
(The number inside parentheses gives the statistical un- 
certainty in the last digits.) A more sophisticated nonlin- 
ear extrapolation of the same lattice data using a chiral 
extrapolation (xE) formula proposed by Leinweber, Lu 
and Thomas (LLT) Q has recently been made by Lein- 
weber and Thomas (LT) [§. The result is -0.16(18). 

This LT estimate contains an unresolved problem 
caused by the use of an extrapolated valence u-quark 
contribution m^o for the baryon "E? that disagrees with 
experiment to renormalize another extrapolated term u„. 
The purpose of this Comment is to suggest that the lat- 
tice data should be renormalized before the xE, and to 
point out that G|/(0) can be extracted more readily from 
the nucleon isoscalar lattice data. The result is -0.55(37), 
where the error comes from uncertainties in the extrap- 
olation formula, in the least-square fitting, and in the 
SU(6) constants F and B. 

To extrapolate a MM term [i(m) calculated on a lattice 
for several large values of the pion mass m to its physical 
value at ttItt, LT use a Pade approximant given by their 
Eq. (13) (denoted here LT13) that is the inverse of a 
quadratic Taylor approximant in m. The extrapolated 
results, needed in this Comment, are read from Figs. 4 
and 5 of 1^ and given in Table |. They agree well with 
experimental values |^ where known, as pointed out in 
for /ip = p and /i„ = n. 

For technical reasons, LT use a simple quadratic Taylor 
approximant in m, Eq. (14) in LT and called here LT14, 
for Mgo . (In terms of s = m} , LT14 is the often- used linear 
approximant improved by a certain leading nonanalytic 
(LNA) term x™--) Table || shows that the extrapolated 
value for it^o differs significantly from the experimental 



value 2*^ — 2^. A concise review of these results in greater 
details has been given by Q. 

The discrepancy from experiment is not due to the use 
of LT14, as I have explicitly obtained an LT13 fit giving 
the same extrapolated value of m^o = -0.37 (and -0.38 
for LT14). The discrepancy comes instead from the sign 
of X that controls the curvature direction caused by the 
LNA term. In the LT model, x has opposite signs for u„ 
and u^o, because the dominant pion cloud has opposite 
signs, namely n — > p7r~ — > n compared to S° — > S^tt"*" 
^ S° Q. This sign makes it difficult for the xE of u~a 
to curve towards the experimental point. This is the 
unresolved problem in LT mentioned at the beginning of 
this Comment 

To extract (dropping the argument for nota- 
tional simplicity), LT use a certain "ratio" method of 
Leinweber under which the disconnected-loop (DL) 
contribution of the nucleon sea is calculated from 
one of the following formulas: 



l^ = oiP + 2n - /„u„). 
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That is, the extrapolated value of u„ is not used directly, 
but only after renormalization by /„. Whenever w^o ex- 
trapolates correctly to its experimental value — S~, 
the correction factor /„ is 1 and therefore has no effect. 
The correction factor differs from 1 only when the ex- 
trapolation is unsuccessful. Then the argument for using 
the nontrivial correction factor /„ ^ 1 is that the ratio 
Un/u^o that appears has smaller systematic errors 

The large renormalization /„ = 1.62 found for Eq. (|l|) 
is the direct consequence of the aformentioned failure of 
the LT extrapolation for u^o. In contrast, fp — 0.973 
needed in Eq. (||) contains little renormalization because 
LT13 works well there. The LT renormalized values (ac- 
tually fiUi) are shown in Table || in the column marked 
LTR with errors from LT. 
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Is this postulated cancellation of systematic errors af- 
ter LT extrapolations real? On the lattice, u„ and u^o 
have very similar values, and their ratios are likely to 
contain smaller systematic errors. Unfortunately, these 
ratios cannot be used directly in the xEs, and the u's 
have to be extrapolated separately for the lattice data 
used. Hence the correct procedure is to renormalize one 
set of lattice data to fit an available experimental MM 
(thus ensuring that fi = l always), and to use the same 
renormalization on the second set of lattice data before 
the xE. 

This scheme can readily be executed for u-^+ since 
the original extrapolation is already close to experiment. 
With the lattice data for different m actually corre- 
lated, moving up and down together, one can simply take 
{ui + xAui} as the corrected input with the unknown x 
chosen to be -0.23 to make the extrapolated value of u-^+ 
tit S+ — E^. The same x is then used to renormalize the 
lattice data for Up. The LT13 extrapolated result, Up — 
4.16, agrees with that from Eq. (||). 

When the same procedure is applied to u-^o , it reaches 
a minimum of only -0.52, a little short of the experimen- 
tal value. (The solution for the input y{ui} scaled by an 
adjustable factor y reaches down to only -0.47.) There 
are LT14 solutions {up = 2.34, u„ = -1.54) that are very 
different from those of the Leinweber equations. Since 
LT14 also does very poorly for the nucleon MMs them- 
selves, giving p = 3.14 and n = -2.39, it appears to be 
unreliable. 

The lattice data can also be renormalized by tilting 
and bending. This is conveniently realized by includ- 
ing the experimental point in an LT13 fit and using the 
resulting set of ratios {ui(output)/Mi (input)} to renor- 
malize the second set of lattice data before its xE. The 
resulting renormalized lattice data (RLD) extrapolate to 
u„ = -0.64(15), showing only a quarter of the renormal- 
ization effect from Eq. (|l|), while Up = 4.13(24) agrees 
well with Eq. (||). These results suggest that Eq. (||) 
is reliable, but Eq. (0) greatly over-estimates the renor- 
malization. The statistical errors shown are obtained by 
treating the lattice data as correlated, i.e., calculated 
from {ui ±pAui}, where p = ^1 - {Np/Nd) = l/\/3, 
where Nd = 3 is the number of lattice data and Np — 2 
is the number of fitted parameters . 

An alternative explanation of the discrepancy in u^o 
that cannot be excluded at the present time is that its 
X parameter is different from that calculated in LT. For 
example, the x parameter needed to fit the experimental 
point together with the original lattice data is x('"h°) 
— 0.2 (-0.3/-I-0.5), showing the change of sign needed for 
the xE to curve towards the experimental point. In other 
words, a relatively small change in the LT value of xi'fJ-E") 
= -0.4 is enough to change its sign. In the rest of this 
Comment, I shall work only with the LT model. 

The value of G%i = -0.16(18) reported by LT is ob- 
tained from Eq.(|l|) using the ratio Un/u^o — 1.51(37). 



(The reported statistical error of 0.18 for G%j is a mis- 
print. It should read 0.27.) LT has also obtained 
C'li = -0.57(42) from Eq.(||) using the ratio Up/u^+ = 
1.14(11), unchanged from the Leinweber linearly extrap- 
olated value of 1.14(8). This is the better estimate be- 
cause the renormalization involved is so much smaller 
and better established. 

The second issue raised in this Comment is that the 
nucleon isoscalar quantity G%j can be extracted from nu- 
cleon isoscalar lattice data alone without using any infor- 
mation from u^o or u-^+ . The resulting physical picture 
for G%j is sufficiently simple to permit a rather unique 
answer to be obtained for any extrapolation formula used 
such as LT13. Rough estimates of other uncertainties can 
then be made readily. 

The quantity to be extracted is the isoscalar DL con- 
tribution that, under isospin symmetry, is simply 

fJ-N^ = '^{P + n)-^{Up + Un). (3) 

The numerical answer is = -0.18(5) using the ex- 
perimental value oi p + n and the unrenormalized LT 
results shown in Table ||. The answer is -0.10(7) with the 
original LT renormalization, and -0.14(9) using RLD to 
renormalize Up and u„. The statistical errors shown are 
obtained from those in Up and u„ treated as independent. 

To do better in both value and error, it is necessary 
to extrapolate Up + u„ together. Fortunately, there is 
enough information in the lattice data tabulated in LWD 
||l| to extract the covariance ((AMp)(Au„)) and to calcu- 
late the statistical error for any linear combination aui + 
hu2 on the lattice. The combination is then extrapolated 
to m-r; by using the chiral parameter axi + bx2- 

The least-square fitting program and input lattice data 
used here are first validated by checking against pub- 
lished results. The top panel of Table || gives a compar- 
ison with the results of LLT Q for the nucleon MMs p 
and n. The data used come from two independent sources 
||l|,|| . My errors are obtained by assuming that all lattice 
data are correlated (i.e., using Np — 2 and Nd — tS). The 
chiral parameters x ^re defined in terms of the SU(6) 
constants F and D ||l^,P . The one- loop corrected values 
|]lO| are used in both LLT and LT. The table shows good 
agreement in both extrapolated values and fitting errors. 

In the second panel of the table, extrapolations are 
obtained for the correlated lattice data |l| used in LT. 
Both values and errors agree with their results. 

The third panel of Table || shows the best results from 
Tabcl ^ obtained by extrapolating Up ± u„ together, 
using different x parameters for flavor and valence con- 
tributions. The quality of this single-step extrapolation 
is first checked for the nucleon MMs. The isovector mo- 
ment p — n = Up - Un on the lattice because the isoscalar 
sea contributes nothing. Table |l| shows that the p — n 
result of 4.65 for LT13 agrees well with the value of 4.68 
from separate extrapolations. Both differ somewhat from 
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the value Up - u„ = 4.81 from separate extrapolations. 
The isoscalar moment extrapolates to p + n = 1.03, in 
agreement with the value of 1.10 from separate extrapo- 
lations and with experiment. In contrast, LT14 appears 
inadequate for both nucleon MMs. 

The much more successful formula LT13 proposed by 
LLT has the correct l/niq oc 1/m^ behavior for heavy 
quarks. Its comparative success in chiral extrapolations 
comes from a reduction of the effect of the LNA term 
when placed in the denominator. Finally, the importance 
of having the correct heavy-quark limit is illustrated by 
using another inverse Taylor approximant M4D where 
the quadratic term in m in the denominator is replaced 
by an incorrect quartic term. It works well too for the 
large isovector moment {p — n)/2 that dominates the nu- 
cleon MMs. This suggests once again that the LNA term 
should be in the denominator. 

The extrapolated value of Up + Un is also needed to cal- 
culate from Eq. (^. It is obtained from the lattice 



data for 



Un by extrapolating with a valence param- 



eter x" different from the flavor value X'^ — used for 
p + n. The results obtained with different extrapolation 
formulas are shown in Table |l|. The very small statisti- 
cal error found for /x^^, denoted in the table as fi^^{l), 
comes from the fact that while each term in Eq. ^ 
varies substantially as all lattice data move up and down 
together, the difference between the two terms remains 
essentially unchanged. 

To determine if p -I- n should be fitted in extracting 
fj.^^, I renormalize the lattice input data to {ui + xAui} 
where x is adjusted to reproduce the experimental value. 
The changes for all three formulas, also given in the table, 
are all quite small, showing that the procedure is quite 
robust. 

Unfortunately, there is a problem in the extrapolations 
for p + n, and also for p and n separately. The lattice 
data used do not contain any DL contribution, which has 
in fact been assumed negligibly small in all LT extrap- 
olations. Unpublished lattice data from Dong, Liu and 
WiUiams(DLW) |l^jl^ seem to indicate that the DL con- 
tribution is not necessarily negligible in the mass range 
771 « 0.6-0.9 GeV of the LWD data. Unfortunately, the 
covariances for these DL lattice data have not been saved, 
thus precluding a quantitative estimate of the missing DL 
contribution for different 777's. This missing DL term is 
formally twice as strong and fractionally 4-6 times more 
important in p + 77 as it is in each MM separately. It 
gives an unknown systematic error to the extrapolated p 

+ 77. 

Given this unresolved problem in the lattice data, one 
might want to use the experimental value of p + 77 in Eq. 
(^. The results, shown in the table as fj,^^{2), can differ 
significantly from {1). The large error shown comes 
from the statistical fitting error of Up + Un- 

These two estimates of fj}^ agree if the extrapolated 
Up + Un agrees with experiment. It is therefore interest- 



ing to see what happens when the lattice data for 7730 
and 77s+ are used to renormalize the lattice data for 7i„ 
and Up before extrapolating 77^ ± Un by LT13 and M4D. 
The results, given in the lower panel of Table |l[ show 
that the extrapolated p + n hy LT13 is now close to the 
experimental value. I take fi^i^.) as the better estimate 
because it contains no systematic error from p + n and 
a more conservative error estimate. This is the "best" 
result shown in Table |. The difference of 0.07 between 
LT13 and M4D shall be taken to be an additional uncer- 
tainty arising from incomplete knowledge of the correct 
XE formula. This large uncertainty has not been included 
in LT. 

Errors also arise from the choice of the SU(6) constants 
F and D on which the xs depend. The x parameter for 
the isoscalar MM p + n is always X'^ — ^ for any F and 
D. The valence parameter for Up + Un is proportional 
to (and has the same sign as) 



= -a{F + Df, 



5 — 6r + 9r 
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(l + r)2 



(4) 



where r = F/D. The factor a is always positive, so that 
both and G%[ extracted from it are always negative 
in LT extrapolations. The factor a has a minimum of 2.39 
at r = 0.74, and the values 2.40 and 2.45, respectively, at 
the one-loop corrected value of r = 0.66 and the tree-level 
value of r = 0.58 ||l0|] . Hence the error from uncertainties 
in r is small and will be neglected. 

With r — 0.66, the uncertainty m F + D can be deter- 
mined by treating it as a free parameter in LT13 fitted to 
both experimental and lattice data for each nucleon MM 
p or n. The expanded data set used by LLT then gives 



(F + i:>)p = 0.96(20), (F + D) 

A(77p + 77„) = 0.21, A/1 



DL 

N 



1.00(20), 
0.04. 



(5) 



This yields a total error of 0.10 for reported in Table 
^ when quadratically combined with the errors of 0.06 
(from least-square fitting) and 0.07 (from the xE for- 
mula) . 

The resulting values of G%j are shown in the lowest 
panel of Table | for iij = 0.55 and 0.59 as used by LT 
and LLT, respectively. The value of -0.55(37) at = 
0.55 agrees with the result -0.36(20) of DLW ^ ob- 
tained from a direct evaluation of the DL contribution. 
The value — 0.55 used here is actually their ratio of 
directly evaluated strange to u/ d DL contributions. 

From the perspective of Eq. (|^), the Leinweber Eq. 
(|l|) or (H) contains a spurious nucleon isovector term 



A/7 



DL 

N 



n - fi{Up - Un)], 



(6) 



where fi — fp or /„. With fp close to 1, this term is very 
small in Eq. (|^). This is another reason why Eq. (^) is 
better than Eq. M). 
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Additional systematic errors not taken into consider- 
ation here include effects from finite-volume and finite- 
spacing errors and from the quenched approximation not 
already simulated by the xE formula used Q . They also 
include uncertainties in the ratio i?^ and differences be- 
tween the LT model and true QCD, such as the contri- 
butions of neglected kaon loops and of m^lnrnq terms, 
in as far as they affect the parameter analogous to /3" of 
Eq. (Q) that drives the DL contribution in the xE. 

I would like to thank Drs. D. Leinweber and K.F. Liu 
for helpful comments. 
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TABLE I. Comparison of extrapolated magnetic moment 
fii (in units of hm) obtained by Leinweber, Lu and Thomas 
(LT) and by Leinweber and Thomas (LT) ^ without and 
with LT renormalization (LTR) with experimental moments 
and with results obtained here using nucleon lattice data only 
(third panel of the table). The number inside parentheses 
gives the statistical uncertainty in the last digits. Estimates of 
GliiRd) in the lowest panel are obtained from defined 
by Eq. (Bl) using the experimental value of p -I- n. 
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-0.55(37) 
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-0.76 
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-0.42 





TABLE IL Extrapolated nucleon magnetic moments ob- 
tained by using the chiral extrapolation formulas LT14, LT13 
and M4D from nucleon lattice data without and with data 
renormalization to fit p -I- n. The lower panel gives results 
obtained from renormalized lattice data (RLD) for Up and u„ 
that use information from Ugo and Uj,+ . 
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M4D 


4.43(26) 


0.70(5) 


3.24(16) 
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0.88 
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RLD: 
LT13 
M4D 


4.67(42) 
4.44(26) 


0.91(10) 
0.67(4) 


3.53(33) 
3.12(15) 


-0.134(4) 
-0.186(4) 


-0.15(6) 
-0.08(3) 



Up + Un 3.53(33) 
p + n 0.91(10) 0.8798... 
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